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ABSTRACT 

Near-infrared observations of stellar orbits at the Galactic Center provide conclusive 
evidence for a massive black hole associated with the compact radio source Sgr A*. 
The astrometric reference frame for these observations is tied to a set of red giant 
stars, which are also detectable at radio wavelengths through SiO maser emission 
in their envelopes. We have improved the precision and long-term stability of this 
reference frame, in which Sgr A* is localized to within a factor 5 better than previously: 
^ 0.17 mas in position (in 2009) and ~ 0.07 mas yr“^ in velocity. This improvement 
is the result of modeling and correcting optical distortion in the VLT/NACO imager 
to a sub-mas level and including new infrared and radio measurements, which now 
both span more than a decade in time. A further improvement will follow future 
observations and facilitate the detection of relativistic orbital effects. 
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1 INTRODUCTION 

The observability of the Galactic Center in the near-infrared 
and its relative proximity at a distance of ~ 8 kpc allow 
the study of a galactic nucleus in unparalleled detail (for 
a review see Genzel et al. 20101. Adaptive optics assisted 


imaging observations of the Galactic Center with NACO 
at the VLT or NIRC2 at the Keck observatory routinely 
achieve a resolution close to the diffraction limit and a 
depth limited by source confusion. Since the discovery of 


stars with large proper motions (Eckart & Genzel 1996 


Ghez et al. 19981, the continued monitoring of individual 


stellar orbits has produced strong evidence for a compact 
mass of ~ 4 • 10®M© that can be convincingly attributed to a 
quiescent massive black hole (e.g. Ghez et al.|2008 Gillessen 
let al.|[2009|| at the location of the radio source Sgr A* 
i MentG et al.|1997 Reid et ah]|2003 20071. Increasing the 
astrometric precision would permit measuring orbits of more 
stars than presently possible, permit refined estimates of the 
black hole’s mass and distance, and eventually permit the 
detection of relativistic effects on the orbits of short-period 
stars, most immediately post-Newtonian effects on the orbit 
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of the star S2 (e.g. 

Jaroszynski 1998 

Fragile & Mathews 

2000| Rubilar & Eckart||20011 jWeinberg et al.||20051 

Zucker 

et al.|2006| Nucita et al. [20071 |Will|2008| Merritt et a 

.120101 

Angelil et al.|2010| Eisenhauer et al. 12011 1. 


While the motions of the so-called S-stars in the 
central 1 " are measured from images taken with the S13 
camera of NACO (~ 13 mas px“^, see Rousset et al.||isi^ 
Lenzen et al.|19^l, less frequent observations with the S27 


camera (~ 27 mas px are used to efficiently set up the 


astrometric reference frame (Trippe et al. 20081. However, 


the S27 camera is affected by optical image distortion on 
a non-negligible level, as for example is the narrow field 


camera of NIRC2 (Yelda et al. 20101. By systematically 


altering the relative positions of stars on the detector, this 
effect becomes a substantial source of uncertainty when 
positions measured on widely dithered images are combined 
into a common coordinate system or when they are aligned 
with astrometric reference sources. 

In our earlier studies of stellar orbits at the Galactic 


Center (Schodel et al. 2002 2003 Eisenhauer et al. 2003 


|2005[ |Gillessen et al.||2009[ ) we used high-order coordinate 
transformations to mitigate the effect of image distortion 
and relied primarily on a reference frame based on the 
assumption that a sufficiently large number of evolved stars 
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2 Plewa et al. 


in the nuclear cluster show no net motion. The stellar 
motions in this ‘cluster rest frame’ were then placed in 
an astrometric coordinate system centred on Sgr A* by 
aligning the positions of several SiO maser stars at a 
certain epoch with their positions as predicted from radio 
observations. The comparison of the S-stars’ positions over 
several years was still limited by a systematic uncertainty 
of ~ 2 mas associated with the definition of the coordinate 
system, likely because of residual distortion in the S27 
camera. Yet the typical uncertainty of a single position 
is as low as ~ 0.3 mas, limited by residual distortion in 
the S13 camera, uncertainty in modelling the point spread 


function (PSF) and ultimately source confusion (Fritz et al. 


2010). The long-term stability of the cluster rest frame is 
also fundamentally limited to a j « 0.07 mas yr“^ by the 
intrinsic velocity dispersion of the selected stars in the plane 
of the sky (a ~ 3 mas yr“^; Trippe et al. 2008) and the 


number of available stars (n ~ 2000; Gillessen et al.|2009 l. 

In this paper we construct an alternative astrometric 
reference frame, which has now become more advantageous, 
by relating many well-measured stars to the SiO maser 
stars directly at multiple epochs and then using their 
motion with respect to Sgr A* for reference. The details 
on the construction of this ‘Sgr A* rest frame’ ( |Yelda| 
et al. polot are described in section and the precision of 


localizing (radio-)Sgr A* in the infrared frame has important 
implications for the analysis of stellar orbits (see section]^. 
To improve this precision we have implemented a more 
accurate distortion correction for the S27 camera of NACO. 
In section the apparent distortion is measured at four 
epochs between 2004 and 2012 from images of globular 
clusters in the ESO archive. This is possible by matching 
them to Hubble Space Telescope (HST) astrometry, which 
can be considered distortion-free in comparison. Based on 
the resulting distortion models it is then possible to also 
correct other images taken close in time. A per night 
self-calibration approach is not feasible for our Galactic 
Center observations, because the requirement that many of 
the same stars are placed in many different regions on the 
detector is not met. To start with, we describe our data and 
basic data processing in section 


2 IMAGE PROCESSING 
2.1 Observations 



Figure 1. Example of a PSF grid obtained from a Galactic Center 
image with median Strehl ratio. Each PSF is normalized to unit 
flux. The adaptive optics guide star is located near the top of the 
image. The typical deviation from the mean PSF in a single pixel 
is about 10% to 35%. 


optics system. The dithering scheme is typically a wide 
square pattern, such that the total field of view is 42" x 42" 
and the central 14" x 14" are present in all images. At each 
pointing position either two, four or six images are taken, 
before moving on. A small offset of a few pixels is applied 
on revisiting a pointing position to minimize the effect of 
detector artefacts on mosaic images. 

We also make use of high-quality HST proper motion 


catalogues available for oj Centauri ( Anderson fc van der 
|Marel|[2010[ ) and 47 Tucanae (Bellini et ah, 2014, private 
communication), which provide us with a nominally 
distortion-free reference system. The HST astrometry profits 
from the stable observing conditions in space and thanks to 
a well-characterized PSF and distortion correction for the 
instrument channels utilized, the measurement uncertainty 
of stellar positions on a single image is as low as ~ 0.01 px 
or ~ 0.5 mas (e.g. Anderson fc King|2006 ). More details on 
the HST observations can be found in Bellini et al. (2014) 
and references therein. 


2.2 Image Reduction 


We make use of two sets of observations with NACO/S27 
for different purposes, which target the Galactic Center, 
but also the densely populated cores of the two globular 
clusters cj Centauri and 47 Tucanae. The globular cluster 
observations that result in viable distortion corrections are 
summarized in table The instrument was configured 
to use the Ks-band filter, the Aladdin 3 detector in 
the high dynamic range mode and the readout mode 
read-reset-read. This is also the typical setup for many 
epochs of Galactic Center observations using the same 
camera, which together establish the astrometric reference 
frame and are summarized in table [2] A small number 
of images are omitted, for example those compromised by 
severe saturation of the relevant SiO maser stars or those 
affected by PSF artefacts caused by a failure in the adaptive 


The standard steps of sky subtraction, division of a flat-field 
and a bad pixel correction are applied to every raw 
NACO image. The nightly sky background is estimated by 
taking the median image of either a set of dedicated sky 
exposures or else the randomly dithered science exposures. 
The flat-field is created by stacking and normalizing 
lamp exposures taken during daytime calibrations. The 
unusable hot or dead pixels are identified by comparing 
the ADU counts in neighbouring pixels, from which a 
replacement value is interpolated if necessary. Individual 
detector integrations were recorded for most images since 
2010, so that frames of exceptionally low quality can be 
rejected before the remaining frames are averaged to create 
the final image. 
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Table 1. Summary of globular cluster observations: A distortion correction for NACO. 



ESO Prog. ID 

Date 

Target 

Filter 

DIT (s) 

NDIT 

Images 

Sample Size 

Indiv. Sources 

(a) 

074.D-0151(A) 

2004-10-12 

47 Tuc. 

Ks 

2.5 

24 

20 

19826 

1698 

(b) 

482.L-0793(A) 

2009-07-03 

(jj Cen. 

Ks 

0.3454 

100 

22 

3910 

534 

(c) 

60.A-9800(J) 

2010-02-07 

oj Cen. 

Ks 

1. 

10 

20 

3658 

388 

(d) 

089.C-0638(A) 

2012-07-22 

47 Tuc. 

Ks 

1. 

30 

10 

9671 

1824 


Table 2. Summary of Galactic Center observations: The Sgr A* rest frame. 


ESO Prog. ID 

Date 

Filter 

DIT (s) 

NDIT 

Pointings 

X Images 

PSF 

FWHM (px) 

Maser 

Stars 

Reference 

Stars 

Distortion 

Correction 

60.A-9026(A) 

2002-04-01 

Ks 

0.5 

8 

10x1 

3.8 

7 

73 

model (a) 

71.B-0077(A) 

2003-05-09 

Ks 

0.5 

120 

19x1 

2.7 

8 

91 

model (a) 

073.B-0085(E) 

2005-05-12 

Ks 

0.5 

60 

16x6 

2.9 

8 

91 

model (a) 

077.B-0014(A) 

2006-04-28 

Ks 

2.0 

28 

8x4 

3.3 

7 

87 

model (a) 

077.B-0014(E) 

2006-08-27 

Ks 

2.0 

28 

8x4 

4.5 

7 

79 

model (a) 

078.B-0136(B) 

2007-03-16 

Ks 

2.0 

28 

7x4 

2.6 

8 

91 

model (b) 

179.B-0261(A) 

2007-03-31 

Ks 

5.0 

6 

16x6 

2.6 

8 

91 

model (b) 

179.B-0261(M) 

2008-04-04 

Ks 

2.0 

28 

8x4 

2.8 

8 

91 

model (b) 

179.B-0261(N) 

2008-08-04 

Ks 

1.0 

57 

8x4 

3.0 

8 

91 

model (b) 

179.B-0261(U) 

2008-09-15 

Ks 

1.0 

57 

8x4 

2.7 

8 

91 

model (b) 

179.B-0261(X) 

2009-03-28 

Ks 

1.0 

60 

8x4 

2.9 

8 

91 

model (b) 

179.B-0261(X) 

2009-03-30 

Ks 

1.0 

60 

8x4 

2.7 

8 

91 

model (b) 

183.B-0100(J) 

2009-09-19 

Ks 

1.0 

60 

16x4 

2.7 

8 

91 

model (b) 

183.B-0100(J) 

2009-09-20 

Ks 

1.0 

60 

8x4 

2.8 

8 

91 

model (b) 

183.B-0100(T) 

2010-05-08 

Ks 

0.5 

126 

8x4 

3.2 

8 

90 

model (c) 

183.B-0100(V) 

2010-09-27 

Ks 

1.0 

126 

8x2 

3.4 

8 

86 

model (c) 

183.B-0100(X) 

2011-04-01 

Ks 

1.0 

66 

8x4 

3.1 

8 

86 

model (c) 

183.B-0100(V) 

2011-05-16 

Ks 

2.0 

9 

8x2 

2.6 

8 

91 

model (c) 

088.B-1038(A) 

2012-03-14 

Ks 

1.0 

30 

41x1 

3.3 

8 

86 

model (d) 

088.B-0308(B) 

2012-05-03 

Ks 

0.9 

33 

8x6 

3.4 

8 

86 

model (d) 

089.B-0162(D) 

2012-08-08 

Ks 

1.0 

60 

8x4 

2.7 

8 

91 

model (d) 

091.B-0081(F) 

2013-05-13 

Ks 

0.9 

33 

16x6 

2.6 

8 

90 

model (d) 


2.3 Star Lists 

With the aim of achieving the highest astrometric precision, 
we extract the detector positions (and fluxes) of stars in 
every reduced NACO image using an empirical model of 
the PSF (see [Anderson King| |2000| ). To identify the 
stars, an image is first correlated with the current PSF 
estimate, with the correlation operator being the normalized 
cosine distanc^ Initially, a Gaussian kernel is substituted 
for the PSF and since this kernel is symmetric, all PSFs 
derived from it are naturally centred. The result is a map of 
coefficients between —1 and 1, which measure at every pixel 
how well the image matches to the PSF locally. The location 
of a candidate star is the peak of a connected pixel region 
in this map, for which the correlation coefficient exceeds 
a threshold value of 0.7 and the ADU counts fall in the 
range between the noise level and the full-well capacity of 
the detector. At each candidate position, the PSF is fit to 
a roughly circular, top-hat image region by method of least 
squares. There are a free position offset and a free scaling 
factor, as well as three additional parameters that describe a 
tilted planar background, which minimizes a potential bias 

^ This distance is . ■ 7~7| between two vectors u and v, 

jli —14| \V — V\ ’ 

which are here one-dimensional arrangements of the correlation 
kernel and the neighbourhood of a certain pixel, respectively. 


in crowded regions. The radius of the fit region, i.e. the 
number of pixels surrounding a central pixel at all sides, 
is typically 4 px, but tied to the full width at half-maximum 
(FWHM) of the PSF such that the entire PSF core is 
included. 

The actual PSF not only changes with time, but 
varies spatially too. For example, the PSF gets significantly 
broader and elongated at distances of 10” to 20” from the 
adaptive optics guide star, as the correction of atmospheric 
turbulence degrades due to anisoplanatism. To take into 
account such a variation we use one 4x4 grid of PSFs per 
image, each of which is derived from up to fifty ideally bright 
and isolated stars within cells of 256 x 256 px size spanning 
the detector (Fig. [^. By using a bicubic interpolation to 
resample the surrounding image regions, the selected stars 
are magnified by a factor 2 and centred in sub-images 
typically 35 px wide, depending on the FWHM of the PSF, 
and also normalized to zero background and unit flux. The 
median superposition of these sub-images, weighted by the 
square root of the original fluxes, is an estimate for the PSF 
at the centre of a grid cell. Afterwards, a modified cosine 
window function is used to taper the extended PSF haloes. 
The PSF at a certain detector position is finally estimated 
by means of a bilinear inter- or extrapolation on the whole 
grid, using spline interpolations of third order to evaluate 
each grid PSF at sub pixel offsets. 
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The two iterative steps of extracting a PSF grid and 
creating a star list are repeated three times. In the case 
of mildly saturated stars, pixels with an ADU count above 
80% the full-well capacity are ignored throughout, as long 
as the fraction of excluded pixels is small (< 10%). Heavily 
saturated stars are never identified as candidate stars in 
the first place. The final outcome are lists which contain 
the positions of all detected stars in the pixel coordinate 
system of each particular image. They also contain the 
fluxes of the stars, which are converted into instrumental 
magnitudes. Only stars brighter than the distribution’s peak 
magnitude are kept (m^ fS 16). These are likely to be 
genuine detections for which positions can be measured with 
similar precision. 


3 A DISTORTION CORRECTION FOR NACO 


3.1 Star Matching 

To measure the distortion in the globular cluster images. 


the respective star lists (Section 2.31 are first matched to a 
proper motion catalogue obtained with the HST. A reference 
star list is created for each observing night by propagating all 
catalogue positions to the right date. Since the images from a 
single night overlap, it is convenient to determine the relative 
offsets between the NACO star lists by cross-correlating the 
original images. Initial transformations that shift and scale 
the HST star list to roughly align with each NACO star list 
are determined manually. Matching pairs of stars are then 
found using a nearest neighbour search. For each detected 
star, to make a match with a catalogue star, there has to 
exist exactly one star within a certain radius around it on the 
transformed HST list and no other star on its own NACO 
star list. The latter requirement effectively excludes stars 
with close neighbours. Additionally, the match needs to be 
unique. Once matching pairs of stars are known, they can be 
used to find better transformations for the original HST star 
list and another matching process follows. The cutoff radius 
for the nearest neighbour search is set to half the FWHM of 
the PSF, but five times larger for the first iteration. A last 
quality criterion is the number of matched stars per image. 
It should be greater than twenty to allow the combination 
of multiple star lists later on. Typically, a few thousand 
detector positions of a few hundred individual stars are 
accepted in total, which cover the detector reasonably well. 


parameters), which already combines translation, rotation, 
scaling and shear operations. The map of transformation 
residuals, i.e. the remaining differences between aligned 
positions, is a vector field that represents the image 
distortion. However, to create a single map from multiple 
images, it is necessary to at least once apply a preliminary 
distortion correction to the NACO star lists, recalculate the 
transformations and refit a model to the pooled residuals. 
This iteration reduces a potential bias that occurs if stellar 
positions from individual lists are not distributed uniformly 
over the detector. 

The distortion model we adopt is a linear combination 
of vector fields from a complete orthonormal basis that can 
describe arbitrary vector fields on the unit circle (Table [AI| . 
Such a basis can be constructed from the union of a 
set of vector polynomials with zero curl and another set 
with zero divergence (see Zhao fc Burge|[2007| 20081. Both 
sets are needed to account for aberrations of the optical 
system and the arrangement of the detector. This model 
is fit to each of the four distortion maps using a least 
squares technique, after the detector coordinates have been 
rescaled accordingly and the maps have been smoothed with 
a median filter to clean them from outliers (e.g. falsely 
matched stars). We pick a smoothing length of 64 px 
to roughly correspond to the scale at which the median 
absolute deviation between the smoothed and the original 
vector fields becomes approximately constant. 

Twenty free parameters are needed to fully capture 
the spatial variability of the image distortion. This number 
can be inferred by inspecting the fit residuals, which 
stop decreasing when enough basis fields of higher order 
have been included. Even so, most of the distortion can 
be attributed to the low-order aberrations y-tilt, defocus, 
coma and spherical aberration. In terms of the number 
of parameters, an equivalent model could be composed of 
two third-order polynomials in two variables, each of which 
describes the distortion along one dimension. 

The characteristic shift required to move a measured 
position into a pixel coordinate system approximately free 
of distortion is ~ 0.2 px (~ 5.4 mas), but the actual shift 
varies strongly with location on the detector and can be 
as high as ~ 0.7 px (~ 19 mas) in the lower left corner. 
The distortion pattern changes distinctly as the pairwise 
rms deviation between subsequent models is on the order of 
~ 0.1 px (~ 2.7 mas), but the general pattern is static. 


3.2 The Distortion Model 

The image distortion we aim to correct is expected to be 
static, but might still change sporadically due to differences 
in the optical alignment before and after interventions 
on the instrument, for example. All our Galactic Center 
observations, spread over the whole time-range of NACO 
operations from 2002 to 2013, are in the end covered by 
four different distortion corrections (Fig. [^& Fig. [^. 

At each of the according four epochs of globular cluster 
observations (Table[^ , a linear transformation is determined 
for every NACO star list that aligns the matched reference 
positions from the HST catalogue with the stellar positions 
as measured on the detector, in a least squares sense. 
More precisely, it is an affine transformation (with six 


4 THE SGR A* REST FRAME 
4.1 Image Registration 

To be able to measure the motions of stars later used as 
astrometric reference sources, individual images from each 
night of Galactic Center observations (Table must first 
be combined into a common coordinate system. 

We start by combining the star lists of consecutive 
images belonging to the same pointing position, which share 
similar observing conditions. Stars are iteratively matched 
between those lists, allowing for small offsets between the 
images (< 0.2 px). New star lists are then created from the 
mean positions, but, to further purge the lists of spurious 
detections, only stars detected in all images are kept. The 
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Figure 2. The evolution of image distortion in the S27 camera of NACO. Shown is the shift required to move a measured pixel position 
into a pixel coordinate system approximately free of distortion (black arrows). Also shown is the number of reference sources in different 
regions on the detector (grey-shaded bins). Each model is only valid for a limited period of time, as shown in Fig. 



year 

Figure 3. Timeline of records in the NACO archive, on which are marked the different epochs of globular cluster observations and 
the periods of time during which the corresponding distortion corrections are valid (in colour), as well as the epochs of Galactic 
Center observations that together establish the astrometric reference frame (in black). The exact dates of significant changes in the 
image distortion might coincide with interventions on the instrument, but cannot be determined exactly from our data. However, it is 
straightforward to decide whether a certain distortion correction is applicable by comparing the resulting image registration accuracy. 


measurement uncertainty of the new positions is estimated 
by the standard error of the mean and has a typical value 
of about ~ 0.3 mas (mode of histogram). 

The single-pointing star lists are then combined into 
one master star list per epoch. AfBne transformations 
should suffice to register the star lists after they have 
been corrected for image distortion, since any remaining 


non-linear displacements should be negligible (considering 
for example atmospheric refraction, aberration due to the 
motion of the Earth, light deflection in the gravitation field 
of the Sun or curvature of the celestial sphere). Nevertheless, 
we explicitly correct the star lists for differential achromatic 
refraction as well, to minimize anisotropy of the pixel scale 
in the final master list (see Gubler fc Tytler|1998 1. 
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Figure 4. Best-fitting models for the measurement uncertainty 
of stellar positions as a function of K-band magnitude (dashed 
lines) and the mean model (solid line), which is applied whenever 
the uncertainty cannot be estimated reliably from the images 
themselves (see text). 


The initial master star list is the first star list containing 
the maximum number of SiO maser stars. We proceed by 
iteratively matching to it the list with the next smallest 
pointing offset, finding an affine transformation between 
the matched positions and creating a star catalogue that 
contains all aligned positions, grouped by star (unmatched 
positions are simply appended). The new master list is 
a collapsed version of this catalogue, i.e. it contains the 
weighted average positions and the mean uncertainties. In 
this way both the catalogue and the master star list are 
updated until all lists are merged. The whole stitching 
process is repeated once, with the only change being that 
every star list is matched to the same intermediate master 
star list. 

A few epochs are treated differently in some ways. For 
observations before 2004, some image reduction parameters 
have to be adjusted for the different characteristics of the 
Aladdin 2 detector, which was replaced during that year. 
For observations with either a random dithering scheme 
or less than three images per pointing, the measurement 
uncertainty of stellar positions cannot be estimated reliably 
from the images themselves. Instead, we apply a mean error 
model based on the other observations by fitting a power 
law to the one-dimensional uncertainties of consistently 
detected stars as a function of flux (Flg.|^. The photometric 
zero-point is calibrated using the 91 primary astrometric 
reference stars (section 4.3 Gillessen et al.|2009l. 


4.2 Testing the Distortion Correction 

After the image registration step, any uncorrected distortion 
becomes apparent in non-random structure of the final 
transformation residuals. This is because the positions of the 
same stars, measured in different locations on the detector, 
should align nevertheless. We can thus check the quality 
of the distortion corrections and identify approximately the 
periods during which a particular solution is valid (Fig. |^. 

Every transformation residual is divided by its expected 
uncertainty, i.e. the combined uncertainty of the two 
respective positions, and sorted into bins of 128 px x 128 px 
size on the detector (Fig. [^. The uncertainties of the 
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Figure 5. The mean (blue arrows) and standard deviation (blue 
circles) of the standardized registration residuals in bins on the 
detector, overlaid for all epochs. A systematic deviation from a 
standard normal distribution (red circles) indicates uncorrected 
distortion or unrealistic uncertainties. 


transformation parameters are negligible, given that a 
few hundred well-measured stars can always be matched 
between two star lists. An additional statistical uncertainty 
of ~ 0.1 mas needs to be added (in quadrature) to each 
stellar position to make the overall distribution of the 
so-standardized residuals approximate a standard normal 
distribution. However, the agreement is best at the centre 
of the detector and worse at the edges, where the remaining 
time-averaged distortion is usually at most ~ 0.2 mas. We 
therefore add a different uncertainty in each bin and then 
average over a star’s detector positions. The distortion 
model is inherently less constrained at the edges of the 
detector, but the actual distortion is also expected to vary 
on a low level even during the night, for example due to 
unstable performance of the adaptive optics system. 


4.3 The Astrometric Reference Stars 


The proper motions of several late-type giant stars have been 
measured relative to Sgr A* directly at radio wavelengths by 


observing their circumstellar SiO maser emission (Menten 


et al. 1997 Reid et al. 2003 20071. The results have 


since been updated with the addition of more data and 
an improved analysis, including in particular a correction 
of differential precession, which leads to an apparent 
rotation (Reid et ah, 2015, private communication). Eight 
of these maser stars lie in the central 7” to 24" and are 
typically inside the NACO field of view. Not included is the 
exceptionally bright star IRS 7, because saturation of the 
detector prevents an accurate measurement of its infrared 
position. As the only red supergiant in the sample, its radio 
position would also have a considerably larger uncertainty 
than the ~ 0.5 mas intrinsic to red giants, because it is 
expected to have a larger SiO maser shell (l?> 4 AU) and 
the variation of the maser emission on time-scales of ~ 1 yr 
is not resolved. 

At each of the 22 epochs of Galactic Genter 
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Figure 6. The infrared positions of the SiO maser stars and the best-fitting proper motion models (in colour). In each panel the top 
plots show the motion relative to Sgr A* and the bottom plots show the motion relative to the radio reference measurements (in grey). 
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■ average difference 

■ IRS g 

■ IRS 10EE 

■ IRS 12N 

■ IRS 15NE 

■ IRS 17 

■ IRS 19NW 

■ IRS 28 

■ SiO 15 


Ax / mas 


AVx / mas yr 




■ coordinate system prior 
M, (our coordinate system) 
^ M, (G09 'combined system') 
^ M, (G09 'maser system') 


X / mas 


Vx / mas yr ' 


Figure 7. Top: Differences between the radio- and infrared-measured linear motions of the SiO maser stars in position (coloured lines 
on the left) and velocity (coloured dots on the right), shown with their Icr-uncertainties (dashed ellipses). The position differences are 
evaluated between 1996 (marked by diamond symbols) and 2016 (unmarked endpoints), with the average reference epoch being 2009.02 
(marked by dots in between). Bottom: The drift motion of the central mass required by the best-fitting orbit of the star S2, using 
our coordinate system (cross-hatched region) compared with the ‘combined system’ and the ‘maser system’ from |Gillessen et al.||2b09[ 
single-hatched regions). The average deviation shown above (in black and grey) indicates how precisely (radio-)Sgr A* is localized in the 
infrared reference frame and is used in the fit as prior information on this drift motion, as indicated here (see section[4.4^. 


observations, the infrared positions of the maser stars taken 
from the master star list (Section 4.11 are aligned with 
the propagated radio positions using a weighted affine 
transformation, thus taking into account the uncertainties 
in both sets of positions. The pixel coordinates (Xpx, i/px) of 
all detected stars are thereby converted into angular offsets 
(*, y) from Sgr A’'j^ The positional uncertainties, including 
both measurement uncertainty and the extra uncertainty 
due to residual distortion (Section 4.21, are propagated 
using a Monte Carlo Bootstrap technique, along with the 
uncertainties of the transformation parameters themselves. 

A total of 10"^ realizations of one transformation are 
generated by applying random displacements to the infrared 
and radio positions of the maser stars, according to their 


^ +x increases to the west and +y to the north. The conversion 
to celestial coordinates is given to good approximation by 
(x,y) = (—Aocos(<5), A(5). 


known uncertainties, while simultaneously resampling the 
pairs of corresponding positions with replacement. The 
artificial transformations are applied to the other stars 
as well, whose pixel positions are repeatedly perturbed in 
the same way. A small number of potential realizations 
are excluded (< 1%), since at least three unique pairs of 
positions are needed to define an affine transformation. 
The uncertainty of each star’s astrometric position is 
finally estimated by the standard deviation of the sample 
of transformed positions, which in turn is estimated by 
the median absolute deviation. This statistic is more 
robust against outliers than for example the rms deviation. 
Likewise, the uncertainty of the detector position of 
Sgr A*, typically ~ 0.6 mas, can be found by inverting the 
transformations and finding the pixel positions that map 
on to the origin. Although the typical uncertainty of the 
maser stars’ detector positions is ~ 0.3 mas, that of the 
astrometric positions is ~ 0.8 mas and evidently dominated 
by uncertainty in the alignment. The main reason is the 
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thinly scattered distribution of the maser stars across the 
field of view, but also the additional uncertainty of the radio 
positions. 

We fit each maser star’s astrometric positions as a 
function of time with both a linear and a quadratic proper 
motion model, separately in x and y (Fig. [^& Table [^. For 
simplicity we use one reference epoch per star instead of one 
per fit and thus correlations between fit parameters are not 
entirely eliminated. The uncertainties of the fit parameters 
are estimated from a Monte Carlo sample created by 
performing fits to the many positions of each maser star 
generated at every epoch, which are additionally resampled 
with replacement in time. The mean reduced of 0-97 
suggests that the positional uncertainties have possibly been 
slightly overestimated. 

Whether a detected acceleration is genuine is decided by 
two criteria. First, the direction of the acceleration must be 
towards Sgr A*, i.e. the negative radial component must be 
statistically significant while the tangential component must 
be insignificant, at the 5a level. Second, the magnitude of the 
acceleration must be smaller than the upper limit imposed 
by the gravitational force the massive black hole associated 
with Sgr A* can exert, assuming that this is the dominant 
force and given that the observed projected separation 
from Sgr A* is the minimum true separation. None of the 
best-fitting quadratic models satisfy both criteria and we 
keep the linear models. 

Another step of calibration is necessary to measure the 
astrometric positions of the S-stars on images taken with the 
S13 camera, because the maser stars are either saturated in 
the deeper exposures or outside the smaller field of view. 
The connection between the two image scales is made by a 
sample of 91 astrometric reference stars, which can always be 
observed together with either the S-stars (using NACO/S13) 
or the maser stars (using NACO/S27). These reference stars 
lie in the central 0.8” to 5” and were chosen to be relatively 
isolated from other known sources (Gillessen et al. 20091. We 
fit proper motion models for the reference stars analogously 
to the maser stars. Since also none of them show significant 
plausible accelerations, we discard the quadratic models. 
Lastly, we inspect the fit residuals separately for individual 
stars and epochs. Particularly in 2013 there is a systematic 
pattern originating from a misalignment of the maser stars, 
but it is consistent with the statistical uncertainties. 

The overall detection rate of both the reference and the 
maser stars is about 97%. Some fainter reference stars were 
sometimes not detected in a pointing when they could have 
been, due to a combination of varying image quality and 
confusion. Non-detections of the maser stars were caused by 
saturation of IRS 17 and IRS 9, as well as the placement of 
IRS 12N close to the edge of the north-east images. 


4.4 The Location of Sgr A* 

The best-fitting proper motion models for the SiO maser 
stars can be compared with the predicted radio motions to 
assess the precision and stability of the Sgr A* rest frame 
over time (Fig. [^ & Table [^. We find that the radio- and 
infrared-measured motions of the maser stars are consistent 
with each other and that the average difference motion in 
units of mas is: 


year 


2006 2008 2010 2012 



Figure 8. Precision and stability of the Sgr A* rest frame as a 
function of time when excluding past infrared observations (a) 
and also either excluding radio observations after 2010 (b) or not 
applying a distortion correction (c). 


Ax{t) RS (-0.16 ± 0.17) -b (-0.04 ± 0.08)(t - 2009.02) 
Ay{t) « (-b0.08 ± 0.17) -b (-0.01 ± 0.07)(t - 2009.02) 


We thus conclude that (radio-)Sgr A*, i.e. the origin of 
the radio coordinate system, is localized in the infrared 
reference frame with a precision of ~ 0.17 mas in position 
(in 2009) and ~ 0.07 mas yr“^ in velocity (~ 2.7 km s“^). 
Moreover, the average velocity differences in radial and 
tangential direction indicate that the infrared reference 
frame shows neither pumping (vr/r) nor rotation (v^/r) 
relative to the radio system to within ~ 7.5 /ras yr“^ as“^ 
and ~ 7.0 /ras yr”^ as“^, respectively. 

To calculate the average of the position and velocity 
differences we have used as weights the uncertainties 
of the infrared quantities, which implicitly contain the 
propagated uncertainties of the radio motions. To estimate 
the uncertainties of the average values we have reused 
the Monte Carlo sample of proper motion fits. The epoch 
2009.02 at which positions are compared is the average fit 
reference epoch, with the weights being the uncertainties of 
the absolute velocities. Because we have also used weighted 
alignment transformations to begin with, the absolute 
average values are expected to be nonzero. 

For comparison, Gillessen et al. (20091 were able to 
localize Sgr A* to within ~ 1.80 mas and ~ 0.33 mas yr“^ 
in their ‘maser system’, applying in particular a distortion 
correction for the S27 camera that only allowed for barrel or 
pincushion distortion and a free optical axis (3 parameters, 
see Trippe et al.||2008|. However, by using a more complex 
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Table 3. Proper motions of the SiO maser stars in the Sgr A* rest frame. 


ID 

to (year) 

a^o (mas) 

yo (mas) 

Vx (mas yr l) 

Vy (mas yr 

IRS 9 

2009.17 

-5676.73 ± 0.23 

-6333.35 ± 0.26 

-2.29 ± 0.11 

1.62 ± 0.11 

IRS lOEE 

2009.00 

-7688.65 ± 0.26 

4186.87 ± 0.25 

-0.53 ± 0.12 

-2.87 ± 0.10 

IRS 12N 

2008.56 

3271.96 ± 0.31 

-6912.22 ± 0.34 

1.88 ± 0.11 

-2.49 ± 0.12 

IRS 15NE 

2008.43 

-1218.31 ± 0.46 

11260.96 ± 0.43 

0.67 ± 0.20 

-5.90 ± 0.16 

IRS 17 

2009.04 

-13145.03 ± 0.46 

5542.86 ± 0.46 

0.44 ± 0.20 

-2.35 ± 0.20 

IRS 19NW 

2009.55 

-14542.15 ± 0.63 

-18483.37 ± 0.75 

0.57 ± 0.36 

-2.38 ± 0.36 

IRS 28 

2009.44 

-10478.15 ± 0.27 

-5853.57 ± 0.26 

-0.86 ± 0.14 

-6.64 ± 0.11 

SiO 15 

2009.28 

12458.47 ± 0.58 

-11060.20 ± 0.61 

3.90 ± 0.23 

2.15 ± 0.24 


Table 4. Alignment of the SiO maser stars (positions are compared at epoch 2009.02). 


ID 

Axo (mas) 

Ayo (mas) 

Avx (mas 

; yr l) 

Avy (mas 

1 yr 1) 

IRS 9 

-0.04 

± 

0.23 

0.44 

± 

0.26 

-0.09 

± 

0.11 

-0.02 

± 

0.11 

IRS lOEE 

0.18 

± 

0.26 

-0.38 

± 

0.25 

0.04 

± 

0.12 

0.05 

± 

0.10 

IRS 12N 

-0.09 

± 

0.31 

0.09 

± 

0.34 

-0.00 

± 

0.11 

0.06 

± 

0.12 

IRS 15NE 

-0.26 

± 

0.48 

0.23 

± 

0.44 

-0.04 

± 

0.20 

-0.00 

± 

0.16 

IRS 17 

-2.37 

± 

0.46 

0.95 

± 

0.46 

-0.32 

± 

0.20 

-0.23 

± 

0.20 

IRS 19NW 

-1.35 

± 

0.66 

0.08 

± 

0.78 

0.32 

± 

0.36 

0.15 

± 

0.36 

IRS 28 

0.21 

± 

0.28 

-0.03 

± 

0.26 

-0.10 

± 

0.14 

-0.07 

± 

0.11 

SiO 15 

0.19 

± 

0.58 

-0.47 

± 

0.62 

0.09 

± 

0.23 

-0.07 

± 

0.24 

Average 

-0.16 

± 

0.17 

0.08 

± 

0.17 

-0.04 

± 

0.08 

-0.01 

± 

0.07 


‘combined system’ that involved a velocity calibration in 


the cluster rest frame (Section IT 

1, the latter value reduced 

to ~ 0.06 mas yr 

Yelda et al. ( 

2010| were able to localize 


Sgr A* to within ~ 0.57 mas and ~ 0.09 mas yr“^ in their 
coordinate system, which is also solely based on maser stars 
and constructed from independent observations with the 
NIRC2 imager at the Keck observatory. 

The largest deviations between infrared and radio 
measurements occur for the maser stars IRS 17, IRS 19NW 
and SiO 15. In the radio, IRS 17 has the least certain proper 
motion and was only observed at four of eight epochs, as was 
SiO 15. In the infrared, IRS 19NW is the faintest maser star 
and its bright neighbour IRS 19 might bias its position. It 
is also the star farthest from the adaptive optics guide star 
IRS 7, with the next closer one being SiO 15. 

The fact that the uncertainty of the infrared location 
of Sgr A* is well dehned can be exploited in the analysis of 
stellar orbits, by setting priors on the allowed drift motion 
of the central mas^ The star S2 is currently the most 
important one for constraining the gravitational potential 
and the distance to the Galactic Center. By following 
in fitting a Keplerian orbit to the 
and spectroscopic data of S2 and 
using our updated proper motion models for the astrometric 
reference stars, we were able to reduce this drift motion 
signihcantly (Fig. [^, without needing a cross-calibration to 
the cluster rest frame. 

^ Another prior at hand is the radial velocity of Sgr A*. 


Gillessen et al. (2009 


available astrometric 


5 DISCUSSION & CONCLUSIONS 


The radio source and massive black hole Sgr A* at the 
Galactic Center can now be placed in the origin of an 
infrared astrometric reference frame with a precision of 
~ 0.17 mas in position (in 2009) and ~ 0.07 mas yr“^ in 
velocity, by aligning the positions of SiO maser stars 
measured at both infrared and radio wavelengths. 

Besides the new data collected in recent years, the 
factor 5 improvement over the comparable ‘maser system’ 


from Gillessen et al. (20091 is mainly the result of a 


better, 20-parameter correction for optical distortion in the 
S27 camera of NACO, which we derived by comparing 
archive images of globular clusters to astrometric reference 
catalogues obtained with the HST. The shifts of stellar 
positions on the detector due to image distortion are on 
the level of a few tenths of a pixel (several mas) and the 
general distortion pattern is complex but static. The origin 
of the observed minor, yet distinct and apparently abrupt 
changes is likely a difference in the optical alignment before 
and after interventions on the instrument. We make our four 
individual distortion corrections publicly availabl^with the 
intention that astrometric studies with a different scientific 
focus will benefit. 

The precision and stability of the Sgr A’*’ rest frame will 
continue to improve steadily, under the condition that future 
infrared and radio observations yield more positions of 
previously observed or newly identihed maser stars (Fig.|^. 
A fundamental limit arises, should the maser emission not 
track the centroid of the photospheric emission, owing to the 
spatial distribution of maser spots in the extended stellar 
envelopes. Nevertheless, the already gained improvement of 


http://www.mpe.mpg.de/ir/gc/distortion 
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the infrared reference frame enables a reanalysis of the S-star 
astrometry, to be presented in a forthcoming paper (Plewa 
et ah, in preparation). 

Precision astrometry at the Galactic Center holds great 
scientihc potential and will eventually lead to the detection 
of relativistic effects on stellar orbits. Given an orbit with 
eccentricity e around a black hole of mass M, at a distance 
Ro, the dominant post-Newtonian effect with an impact on 
astrometry is Schwarzschild precession. This effect causes an 
apparent apocentre shift per revolution of ( Weinberg|1972 1: 

GttG M, 


As 


Ro{l-e) 


An extended mass distribution would cause a counteracting 


Newtonian precession, but is still very uncertain (e.g. Genzel 

_ 

0.05 mas yr ^ 


et 


al.||2010p. In the case of the star S2, a drift of the 

would amount to the 
0.8 mas) after one 


coordinate system of' 
magnitude of the former effect (As ^ 
orbital period. The ability to detect a relativistic precession 
of the orbit of S2 therefore hinges on a very stable reference 
frame, but also depends crucially on the astrometric data 
obtained at and around the time of the next pericentre 
passage in 2018. 
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APPENDIX A: THE DISTORTION MODEL 


Table Al. Explicit form of the distortion model in terms of its 
basis vector fields. For a derivation see|Zhao fc Burg^||2007||2008[|. 


V 

Vx{x,y) 

Vy{x,y) 

S2 

a2 

0 

S3 

0 

03 

S4 

y/2a4X 

72041 / 

S5 

V^asy 

y/2a3X 

Se 

\/2aQX 

-i/2aey 

Sr 

y/Qa^xy 

ij\ar{x'^ + 3y^ - 1 ) 

Ss 

ij~^aa(3x'^ + 1 /^- 1 ) 

\f&aaxy 

S9 

2y/Sagxy 

7309 ( 0 : -y)(x + y) 

Sio 

y/3aio{x -y){x + y) 

- 273010 x 1 / 

Sii 

2aiia:(3a:^ + 3y^ — 2 ) 

2oiii/(3x2 + 3i/2 - 2) 

<S’l2 

2\/2cii2^(2ic^ — 1) 

2720121/(1 - 2i/2) 

Sl3 

2\/2ai3y(3x'^ + y'^ - 1) 
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Sl4 

2ai4{x^ — 3xy^) 

2oi4i/(i/^ - 3x2) 

Sl5 

-2ai5y(y2 _ 3x2) 

2oi5(x 3 - 3x1/2) 

T4 

V2b4y 

— y/2b4X 

Tr 

yj^brix'^ + 3y^ - 1) 

—y/Qbjxy 

Ta 

VSbaxy 

-^§68(3x2 + 7-1) 

Til 

2biiy{3x^ + 3j/2 - 2) 

—26iix(3x 2 + 3t/2 — 2 ) 

Ti2 

2726121/(1 - 2i/2) 

2726i2x( 1 - 2x2) 

Ti3 
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